Density-dependence of functional development in spiking cortical networks grown in vitro 



Michael I. Ham 

Center for Nonlinear Studies and Applied Modem Physics P-21, 
Los Alamos National Laboratory, Los Alamos NM 87544 Tel.: +1-505-565-6230 mikeh@lanl.gov 

Vadas Gintautas and Marko A. Rodriguez 

Center for Nonlinear Studies and Applied Mathematics and Plasma Physics, 
Los Alamos National Laboratory, Los Alamos NM 87544 



oo 
o 
o 

(N 
> 

o 

(N 

u 

d 

•i-H 

cr: 



> 

00 

m 



oo 
O 



X 



Ryan A. Bennett 

Center for Network Neuroscience and Department of Physics, University of North Texas, Denton TX 76203 

Cara L. Santa Maria 

Center for Network Neuroscience and Department of Biology, University of North Texas, Denton TX 76203 

Luis M.A. Bettencourt 

Theoretical Division, Los Alamos National Laboratory, 
Los Alamos NM 87544 and Santa Fe Institute, 1399 Hyde Park Road Santa Fe NM 87501 

During development, the mammalian brain differentiates into specialized regions with distinct functional abil- 
ities. While many factors contribute to functional specialization, we explore the effect of neuronal density on 
the development of neuronal interactions in vitro. Two types of cortical networks, dense and sparse, with 50,000 
and 12,000 total cells respectively, are studied. Activation graphs that represent pairwise neuronal interactions 
are constructed using a competitive first response model. These graphs reveal that, during development in vitro, 
dense networks form activation connections earlier than sparse networks. Link entropy analysis of dense net- 
work activation graphs suggests that the majority of connections between electrodes are reciprocal in nature. 
Information theoretic measures reveal that early functional information interactions (among 3 cells) are syner- 
getic in both dense and sparse networks. However, during later stages of development, previously synergetic 
relationships become primarily redundant in dense, but not in sparse networks. Large link entropy values in the 
activation graph are related to the domination of redundant ensembles in late stages of development in dense 
networks. Results demonstrate differences between dense and sparse networks in terms of informational groups, 
pairwise relationships, and activation graphs. These differences suggest that variations in cell density may result 
in different functional specialization of nervous system tissue in vivo. 



INTRODUCTION 

The mammalian brain is a remarkable structure composed 
of many specialized regions and types of cells. Despite or- 
ganizational differences in neural tissue, the basic functional 
units of the nervous system, neurons, are generally simi- 
lar across tissues, as are methods of forming and modifying 
synaptic connections between them (e.g. spike-timing depen- 
dent plasticity ifisll . Thus, functional specialization of 
brain regions is a function of neuron specialization, (e.g. ex- 
citatory or inhibitory), ratios of neurons to neuroglia, synaptic 
density, learning, and many other factors that are not yet well 
understood. 

Though the functional role of neurons in different types 
of tissues can be similar, neuronal density can vary dramat- 
ically. For example, neuronal density in human fascia dentata 
is w 3.2X10^ neurons/mm^ |16iil while cortical tissue density 
is w 3.4X10^ neurons/mm^ ||2|], nearly an order of magnitude 
difference. 

A previous study of network development in vitro demon- 
strated that network bursting (when a large majority of neu- 
rons fire in a coordinated pattern) and spiking patterns are af- 
fected by neuronal density Ii28il . Since spike and burst activity 



are observable results of neuronal interaction (functional con- 
nectivity) they can be used to infer neuronal relationships 1,2 ij] , 

fi], il, 11- 

In this work the effect of density on functional units is 
explored using dissociated cortical tissue developing in vitro 
on microelectrode arrays. These are well established models 
of neuronal interaction ll22ll . ifTHl . though obvious con- 
straints and limitations must be considered when attempting 
to extrapolate between in vitro and in vivo structures iQl, lfl9ll . 
To gain new insight into functional connectivity in developing 
networks, we analyze the coordinated electrophysiological ac- 
tivity of groups of two and three spike trains, each represent- 
ing an integration of all of the action potentials recorded at a 
single electrode. Each electrode may capture the activity of 
a single neuron, or less frequently, incorporate signals from 
several cells. 

Activation between pairs of electrodes is infeiTed using a 
competitive first response model, whereby directional links 
are derived from spike train data. These pairwise links esti- 
mate the probability that activity at one electrode causes ac- 
tivity at another. The set of all such links forms a network 
activation graph where vertices represent recording electrodes 
and weighted edges represent a dependent activation probabil- 
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ity. We use entropy based link analysis (link entropy) to char- 
acterize connectivity between electrodes. Information theo- 
retic measures applied to ensembles of 3 electrodes reveal 
functional information structures. Such structures are char- 
acterized using an information overlap method that reveals 
whether interactions are synergetic (more information is ob- 
tained from the measurement of two electrodes together con- 
ditioned on the third compared to measuring them separately), 
redundant (less information is obtained from the measurement 
of two electrodes together conditioned on the third compared 
to measuring them separately) or independent. Results show 
that functional structures in networks are strongly influenced 
by neuronal density and suggest that varying cell density is a 
potential strategy for differentiating tissue functionality. 



Network Electrodes Recording days 



Dl 56 4-7,9-26,28,31-35 

D2 56 4-26,28,31-35 

D3 56 4-26,28,31-35 

D4 56 4-26,28,31-35 

51 58 4-7,10,11,13,14,17,19,21,22,24-26,28,31,33,34 

52 58 4-7,10,11,13,14,17,19,21,22,24-26,28,31,33,34 

53 58 4-7,10,11,13,14,17,19,21,22,24-26,28,31,33,34 

54 58 4-8,10,12-14,17,18,20,21,31-35 



TABLE I: Network: D1-D4 are densely seeded and S1-S4 are 
sparsely seeded. Electrodes: number of electrodes with action poten- 
tial activity. Recording days: days in vitro a network was recorded. 



METHODS 

We analyze developmental activity patterns in eight cul- 
tured cortical networks growing on microelectrode arrays. 
Networks were recorded t28i1 on days ranging from 3 to 36 
days in vitro (DIV, days after plating). Half the networks stud- 
ied are sparsely seeded (625 cells/ml) the other half densely 
seeded (2,500 cells/ml). 

Cells (neurons and support glia) are removed from embry- 
onic mice and their existing tissue structure is dissociated 
mechanically and enzymatically. After cells are seeded on 
a microelectrode array, they form new connections and self- 
organize into spontaneously active networks lll2|] . OlTll . Figure 
[U depicts a representative dense and sparse network at 1, 15 
and 32 DIV. Experiments and data collection were performed 
by Wagenaar et al. 128(] and made publicly available for analy- 
sis. Detailed information about culturing, plating, and feeding 
techniques are provided in this reference. 



data presented here, typically no more than 2 neurons are ob- 
served at any given electrode 11 1311 . 

The first response model provides a method for estimating 
whether a spike produced at one electrode initiates a spike at 
another Each time a spike at an electrode is the first to fire 
within 1 — 10 ms after a spike at another electrode (Fig. |2^), 
the link (edge) from the first to the second is incrementally in- 
creased. This results in a weighted and directed edge between 
electrodes (Fig. |2j5). 

Edge weights are calculated for all electrode pairs. Edges 
represent inferred activation links in the network and are not 
necessarily representative of actual synaptic connection be- 
tween neurons recorded at any two electrodes. Edges are de- 
noted by Xij , where / is the first electrode to fire and j the 
second. When pairwise edges are studied individually, they 
are normalized by the lesser of the total spikes recorded at / 
or j, giving an estimate of the probability that / leads j. Nor- 
malized frequency edges Xij and xji are denoted X„ and Y„ 
respectively. 



First response model 

We build a model of network connectivity using spike tim- 
ing correlations to study the evolution of pairwise neural in- 
teractions during development. The model is based on the 
assumption that, within a small network, all activity preced- 
ing a spike within a biologically plausible time window con- 
tributes to its firing. Any given spike is assumed to be corre- 
lated with the ignition of the next spike to fire within 1 — 10 
ms. This temporal window is selected such that the first spike 
has enough time to influence the production of the second, 
but not so long that its effect will have faded [20], 01 • All 
spikes collected at a single electrode form a single spike train. 
Therefore, the number of spike trains is equal to the number of 
active channels (Table|l]i and functional connections discussed 
here are made between electrodes. 

Recent analysis of similar recordings in which individ- 
ual neurons at an electrode were discriminated (spike sorted) 
shows that, for plating densities comparable to those in the 



Activation graph 

The combination of all weighted directed edges between 
electrodes forms an activation graph, which is a weighted di- 
rected graph, with each edge normalized between zero and 
unity. This graph estimates the activation probability between 
all electrodes. New activation graphs are created for each net- 
work on each recording day. 

Information ensembles of order 2 

Neuronal activation along pathways in the brain is generally 
a sequential process. Activity in one area can i gnit e activity 
in another on large scales (tissue structure, e.g. 12611 ) as well 
as small scales (neuronal networks, e.g. [13|]). To characterize 
the influence a target electrode has on all other electrodes, we 
apply a link entropy analysis of the activation graphs. Link 
entropy is calculated for each electrode and measures the un- 
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FIG. 1: Neurons and neuroglia growing on microelectrode arrays. Dense (left column) and sparse (right column) networks are depicted at 
1, 15 and 32 days in vitro (DIV). Dark circles are electrodes and dark lines are integrated wires that transmit action potential activity from 
electrodes to amplifiers (not pictured). Figure was previously published by Wagenaar et al. 1280 ^nd is used with permission from the authors 
under the license of the publication 



certainty of which electrodes will be activated by the target 
electrode. For example, in Fig. |2j), electrode X activates only 
Z. Since there is no uncertainty about which electrode X ac- 
tivates, the link entropy is 0. If X activated more than one 
other electrode, uncertainty in activation and link entropy in- 
crease. In a fully connected network with electrodes, there 
are — 1 potential pathways (edges) originating from each 
electrode. Edges (jc,^) originating from a target electrode (/) 
in the weighted directed activation graph are used to estimate 



the probability p{xij): 

Pixij) - 

The link entropy of the target electrode (;) in an activation 
graph with N active electrodes is: 

N-l 

^/ = - L Pixij)l0g2{p{xij)) (2) 
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< 1 ms. 



FIG. 2: a) Activation links are made when activity at an electrode is 
the first to occur within a 1 — 10 ms window after another electrode. 
Note that no link is made from the second spike in X since the time 
window between it and Z is less than 1 ms and it and Y is greater than 
10 ms. b.) Pairwise activation graph constructed from data in part a. 
Edges are directed and weights increase as connections are observed. 
Edges are normalized by the lesser of the total spikes observed at 
either of the connected electrodes. Activation networks represent 
activation flow and are not representative of the physical connections 
between two electrodes 



theoretic quantities are measured only during network bursts, 
where the vast majority of all coordinated interactions occur. 
More details about network burst detection are given in ifisll . 
The average number of detected network bursts per minute 
across all dense and sparse networks on all recorded DIV is 
shown in Fig. |4^. Bursts are digitized by dividing time into 
10 ms bins. If activity occurs during a bin, the bin is assigned 
a 1 . Otherwise the bin is assigned a 0. Probabilities are com- 
puted by normalizing the number of spiking and not spiking 
bins by the total number of bins. Joint probabilities can also 
be computed in the same way for multiple electrodes. 

After digitization, all groups of spike train triplets are ex- 
amined for redundant or synergetic relationships per the con- 
vention established in (|5|], 124)1 . Specifically, for three elec- 
trodes X, F, Z, we consider the conditional mutual information 
between X and the state of two other electrodes Y and Z (19J): 



a. 
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FIG. 3: a) Synergetic relationship where X and Y together provide 
more information about the state of Z than they do individually, as 
in the case of a logic gate, b) Redundant relationship, where X and 
Y together provide less information about the state of Z than they do 
individually, as in the case of a Markov chain 



This measure represents the uncertainty (in bits) of the po- 
tential activation paths originating from a target electrode. 
Link entropy values range from (no uncertainty; activa- 
tion flows along only one path with probability of 100%) to 
a maximal value of log2{N — 1). The latter occurs only when 
(non-zero) probabilistic connections between the target and 
all other — 1 electrodes in the network are equal. 



Informational ensembles of order 3 

Information theoretic measures applied to ensembles of 3 
electrodes are used to estimate functional relationships, be- 
yond simple directed connectivity, in developing networks. 
At this point, probabilities are no longer based on activation 
graph links between electrodes. Rather, each electrode can 
assume two states - spiking or not spiking. In the conven- 
tion established in Reike et al. [23,1 , see also 15J, information 



p{x)p{y,z) 



(3) 



Subtracting this quantity from the sum of I{X;Y) and 
I{X;Z) reveals the informational nature of the ensemble fl. 
For example, in a synergetic ensemble (see Fig. |3^), more in- 
formation is gained when considering Y and Z together rather 
than separately: I{X\Y,Z) > I{X;Y) +I{X;Z). However, in a 
redundant ensemble (see Fig. |3}'), less information is gained 
when considering Y and Z together rather than separately: 
I{X;Y,Z) < I{X;Y) +I{X;Z). If Y and Z are independent, 
I{X;Y,Z) =I{X;Y) +I{X;Z). Therefore, we can use the mea- 
sure/?(X,y,Z) =/(X;F)+/(X;Z)-/(X;F,Z) to characterize 
the functional informational nature of an ensemble. R can be 
positive (redundant), negative (synergetic) or zero (indepen- 
dent). 

To test for statistical significance, experimentally obtained 
R values were compared to values obtained from a Poisson 
null model with the same spiking rates as in the experiment 
jst]. In this model, experimentally obtained spike trains are 
used to form a randomly distributed Poisson spike train for 
each electrode. This eliminates all correlations between spike 
trains, so that non zero values of R are the result of finite sam- 
pling. The Poisson spike train contains the same number of 
spikes as the experimentally observed spike train. In practice, 
R values obtained from this model are very small since coin- 
cidences are minimal between random Poisson spike trains. 
The Poisson model establishes a background noise level be- 
low which R values were considered to be effectively 0. This 
null model was used to establish the background noise level 
below which R values were considered to be effectively 0. The 
mean Poisson R values over 10 instances of the model were on 
the order 10^^ while typical values from the experimental data 
ranged between —10^^ and 10^'. 
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FIG. 4: Comparison of network bursts and pairwise activity in dense and sparse networlcs. Four dense (diamonds) and four sparse (circles) 
are averaged togetlier. a) Average number of network bursts per minute per day for dense and sparse networks. Multiple network bursts per 
minute occur around 9 DIV in dense networks and 22 DIV in sparse networks. Days with no network bursts are not plotted, b) Example of a 
normalized pairwise interaction. X„ is the probability that neuron Y is the first to respond to neuron X and vice-versa for Y„. For all coordinates 
not equal to [0,0], the coordinate magnitude and angle of deviation from the diagonal X„ = Y„ (solid line) is computed, c) Time evolution of 
the average magnitude across all study networks in dense and sparse networks. Note that average magnitudes are relatively small due to the 
fact that very few pairs lead one another a majority of the time. Magnitude in dense networks increases faster, but sparse networks exhibit 
higher vales after 30 DIV. d) Average angle (degrees) of deviation in dense and sparse networks. Throughout maturation, sparse networks show 
greater deviation from the diagonal than dense ones. In both cases deviation decreases with age, but deviation in dense networks decreases 
fasters 



Results 



We analyze eight developing networks, four dense (2.5 ± 
1.5 cells/mm^) and four sparse (0.6 ±0.24 cells/mm^). For 
each network on each recording day dTJi, new activation graphs 
are generated by the competitive first response model (see 
Methods) and directional activation links between pairs are 
established. Thus, pairs can be represented by the coordinate 
[X„ , Y„], where X„ is the estimated probability that electrode X 
activates Y and Y„ is the estimated probability Y activates X 
(Fig|4j)). If Xn=Yn, the point falls on the line y = x (solid line. 
Fig. Each coordinate has a magnitude (distance from 

zero) and a deviation angle from the line y = x. If X„ = Y„ , 
then the deviation angle = and the mag nitude = ^X^ + Y^ . 
In figures & |4}l we examine the magnitude and deviation 
of all non-zero (X„ 7^ & F„ 7^ 0) coordinates in all networks. 
In figure |4j;, the average coordinate magnitude for dense and 
sparse networks on each recorded DIV is plotted. At early 
ages, the average coordinate magnitude in dense networks is 



greater than in their sparse counterparts. This indicates that 
pairwise connections form faster in dense networks. How- 
ever, in sparse networks, the average coordinate magnitude 
stabilizes at (slightly) higher values than in dense networks 
(Fig. Hj;). In FigHJl, we examine the average deviation from 
the line y — x in dense and sparse networks, for all active 
pairs. Pairwise activation is relatively balanced during dense 
network maturation (X„ « Y„), as demonstrated by low devia- 
tion angles. In sparse networks, deviation angles decrease, but 
remain larger than dense deviation angles throughout matura- 
tion. 



Wagenaar et al. Il28ll observed that sparse networks are 
slower to develop network bursting than their dense counter- 
parts (Fig. 2^). Additionally, activation graphs reveal that 
late onset bursting, a function of pairwise interactions 1 13], is 
highly correlated to the slow development of pairwise activa- 
tion connections. Therefore, early development of pairwise 
interactions in dense networks appears to be related to the on- 
set of multiple network bursts per minute around 10 DIV. Note 
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that multiple bursts per minute are not a present in sparse net- 
works until about 22 DIV (Figgli). 



link entropy 

Figure |5] shows the average daily Unk entropy for all elec- 
trodes in each of the eight networks studied. Note that a 
similar number of active electrodes were observed in all net- 
works (Table nil. As dense networks mature, average individ- 
ual electrode link entropy values approach the maximum pos- 
sible value, which indicates nearly uniform probabilistic links 
between any target electrode and all other active electrodes. 
In other words, given activation at a particular dense network 
electrode, the next electrode to activate is almost completely 
random. Mature sparse networks exhibit midrange link en- 
tropy values, indicating that probabilistic connections from 
a target to all other electrodes are unequally distributed and 
therefore more predictable. This corresponds to the larger de- 
viation angles observed in sparse networks (Fig. |4}l). Note 
that sparse network S2 has larger link-entropy values than SI, 
S3 or S4 though not as large as in the dense networks. 



synergetic after 25 DIV. While three out of the four sparse net- 
works remain primarily synergetic (Fig|6] SI, S3, S4), in one 
network (S2) redundant groups reemerge around 30 DIV and 
are present throughout the remainder of the experiment. Note 
that synergetic triplets still maintain a strong presence during 
this time period. This is in clear contrast to dense networks 
where synergetic groups essentially disappear altogether in fa- 
vor of redundant ensembles. This shift is primarily due to the 
evolution of individual triplet R values changing from syn- 
ergetic to redundant. In Figure [T] several randomly selected 
representative triplets and their R values on different record- 
ing days are shown from network D4. We observe that the vast 
majority of dense network triplets (98%) that are synergetic at 
a young age become redundant during maturation. 

At around 21 DIV, redundant ensembles begin emerging in 
both dense and sparse networks (Fig. |6|l. This correlates to the 
time period when cortical networks in vitro are generally con- 
sidered mature lIlSll and when Wagenaar et al. Ii28i] observed 
that spike rates stopped increasing and leveled off. 



Triplet interactions 



Informational relationships of order 3 

Groups of three electrodes (triplets), are the smallest (and 
most abundant) arrangement capable of providing information 
about functional connectivity in a network 101 • As described 
in Methods, informational structures are determined by com- 
paring the information from two electrodes together condi- 
tional on the third with information gained from the two elec- 
trodes separately conditioned on the third. The resulting value 
(R) categorizes activity between electrodes as redundant, syn- 
ergetic or independent. 

On each recording day, R values for all unique electrode 
triplet ensembles are calculated (see Methods). Note that 
R{X,Y,Z) =R{X,Z,Y). In FigurelH distiibutions of R values 
from each network on each recorded day are shown. Syner- 
getic values {R < 0) are grouped in bin sizes of .001 bits while 
redundant values (R > 0) are grouped in bins of size .01 bits. 
The asymmetry in relative sizes of positive and negative R val- 
ues has been demonstrated previously [5]. Independent values 
(R = 0) account for 35% of all R values and are not shown. As 
described in Methods, all values are obtained during network 
burst events to provide a comparison with earlier studies. 

In all instances, young network bursts (10 — 20 DIV) are 
dominated by synergetic triplets. However, around 18 DIV, all 
networks show a shift towards redundant triplets. By 30 DIV, 
matured dense networks are dominated by redundant triplet 
ensembles. This finding is commensurate with our previous 
study of a mature 45 DIV network |^J. Sparse networks are 
observed to switch back and forth between primarily syner- 
getic and primarily redundant from 1 8 — 25 DIV. Unlike their 
dense counterparts, sparse network triplets became primarily 



Large link-entropy values and the presence of redundant en- 
sembles do not appear to be directly correlated to one another. 
However, large values of link entropy seem to be an indica- 
tor that redundant ensembles will dominate the informational 
relationships. It should be noted that in Dl, the rise in aver- 
age link-entropy and the increase in redundant triplets are well 
correlated, but this is not observed in the other three dense 
networks where link-entropy maximizes before redundant en- 
sembles fully dominate. Regardless, results suggest that high 
values of link entropy are necessary for the emergence and 
dominance of redundant neuronal ensembles. Figure|5]S2 in- 
dicates that, with moderate link entropy values, synergetic and 
redundant groups coexist. Relatively low link entropy, like 
that seen in SI, S3, and S4, appears indicative of synergetic 
relationships (Fig. |5]l. 



DISCUSSION 

In this manuscript, we explored quantitative patterns of 
neuronal interactions in developing dense and sparse cultured 
cortical networks. Such networks are formed when dissoci- 
ated prenatal cortical tissue is plated on microelectrode ar- 
rays and new connections are spontaneously formed. Mul- 
tisite electrophysiological data from these networks provides 
unique access to the development of cortical cultures. It 
should be noted that in vitro network formation is guided by 
similar mechanisms (synaptogenesis; il4ll ) used by intact ner- 
vous systems. 
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FIG. 5: Average daily activation link entropy and standard error for dense (Dl — D4) and sparse networks (51 — 54). Averages are plotted as a 
percent of the maximum link entropy, log2{N — 1). Average values in dense networks come close to maximal value, indicating that all nodes 
in the network have nearly equally probable connections with all other nodes. Larger, but not maximal link entropy in S2 correlates to the 
development of redundant groups in Fig. 6. The values plotted are normalized by the maximum possible value, which corresponds to 1 on the 
y axis 



Pairwise interactions 

At the most basic level, network formation can be thought 
of as a pairwise phenomenon where two neurons form synap- 
tic connections with one another through electrical and chem- 
ical means. We extrapolate functional pairwise connections 
through a competitive first response model which is based on 
the assumption that all previous activity (within a biologically 
plausible time window) contributes to action potential initia- 
tion. This model is used to create directed, weighted activa- 
tion graphs. 

Mature activation graphs of dense networks reveal that con- 
nections between electrodes are mostly reciprocal (X„ = Y„, 
Fig|4}J) as indicated by relatively low deviation angles. In con- 
trast, electrodes in sparse networks tend to develop activation 
connection strengths that are skewed in one direction or the 
other (larger deviation angle, X„ ^ Y„, Fig. |4}i). Additionally, 
link entropy in sparse networks is smaller than link entropy in 
dense networks (S1-S4, Fig. |5]). Together, these results indi- 
cate that electrodes show activation pathway biases (low link 
entropy) in sparse networks, but these biases are not recipro- 
cal between pairs. In dense networks, the opposite appears 
true; probabilistic activation pathways between electrodes are 



nearly uniform (large values of link entropy. Fig. |5]l and re- 
ciprocal between pairs. 

In networks processing information, each type of connec- 
tion structure has advantages and disadvantages. Dense net- 
works appear to have more redundancy (high link entropy) 
which may allow for higher fault tolerance: if a path via a 
given electrode becomes unresponsive to activation attempts 
by another, many other pathways could easily be activated. 
However, such a system may be detrimental to rapid infor- 
mation processing where fewer strong pathways, like those 
seen in sparse networks, could help an animal reach a fast de- 
cision state with greater precision. The presence of favored 
activation pathways in sparse networks may indicate that they 
are better platforms for performing training studies that seek 
to change network interactions in vitro |25|] . Since there 
are fewer probabilistic pathways for modification, changes to 
these connections will likely produce greater differences in 
existing efferent and afferent activation pathways. Conversely, 
shutting down or building pathways in dense networks would 
likely produce a less noticeable change in directional activa- 
tion patterns as many detours are available. 

Ensembles of 3 nodes are the minimum required to form 
computational groups |6]. Sequential chains of redundant 
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FIG. 6: Computed R values for all ensembles of 3 electrodes in developing networks. Independent R values (R = 0) are not depicted. Dense 
networks (D1-D4) are depicted on the top row and sparse (S1-S4) on the bottom row. Individual triplets in dense networks change from 
primarily synergetic relationships to primarily redundant at about 18 DIV and by 30 DIV redundant ensembles dominate. At early and late 
stages of maturation, sparse networks SI, S3 and S4 are dominated by synergetic relationships. Starting around 18 DIV, sparse networks switch 
between primarily redundant and primarily synergetic until about 25 DIV when they settle on primarily redundant. Note in S2, redundant 
ensembles that develop around 28 DIV though synergetic relationships do not disappear as in dense network 




FIG. 7: Daily R values of four representative triplets in network D4. Note that some triplets switch back and forth between synergetic (R < 0) 
and redundant (R > 0) informational states during development. However, after 28 DIV the vast majority of network triplets in D4 were either 
redundant or independent 



cells, shown in Figure can be used to relay information, 
and similar arrangements have been shown to play a role in 
short-term memory in the brain fioll . In the synergetic 
configuration depicted in Fig|2^, a neuron receives inputs from 
many other cells in such a way that spiking activity is a nonlin- 



ear function of the inputs. Both systems (redundant and syn- 
ergetic) allow neurons to integrate activity from many sources 
and process information. 

We demonstrated that plating density has an important ef- 
fect on informational relationships. Dense and sparse net- 
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works develop similar distributions of synergetic informa- 
tional structures between 10 and 20 DIV (Fig. |6l). However, 
as dense networks mature, most triplets change from syner- 
getic to redundant functional informational connections (Fig. 
|7|. By about 30 DIV very few synergetic triplets remain in 
dense networks. Conversely, individual triplets in sparse net- 
works remain primarily synergetic from early to late stages 
of development. These results suggest that sparse networks 
may be better suited for processing information since syner- 
getic relationships are necessary to perform these tasks jjstl- 
Of course synergetic relationships also exist in denser net- 
works, but in a smaller relative proportion. They may also 
occur in larger groups of cells, not analyzed here One 
mature sparse network, S2 (Fig. |6]l, showed a relatively even 
mix of both synergetic and redundant ensembles during late 
stages of maturation. This was not seen in any of the other 
networks. Furthermore, this network exhibited link entropy 
values lying between the other sparse networks and the dense 
networks (Fig|5]l. These findings suggest a link between ac- 
tivation connections and informational relationships. Specifi- 
cally, when activation connections are balanced throughout a 
network (larger link entropy values), redundant ensembles re- 
sult. However, if activation connection are not quite balanced 
(mid range link entropy values), but are close, both synergetic 
and redundant ensembles appear 

The observed density dependent functional organization of 
cortical tissue in vitro raises interesting questions for future 
research. To elucidate how network density affects connectiv- 
ity structures, experiments that span a broader range of plating 
densities are needed. Additionally, it would be interesting to 
apply these analytical techniques to recordings where a net- 
work is monitored continuously during development. 



Conclusion 

The development of functional connectivity in neural net- 
works in vitro comprising similar cells is greatly affected by 
plating density. In a comparison of dense and sparsely pre- 
pared networks, we demonstrate that strong pairwise connec- 
tions occur earlier in dense networks. Link entropy analy- 
sis of pairwise activation graphs reveal that nodal connections 
tend to be biased in sparse and nearly equal in dense net- 
works. Connections between electrode triplets in both dense 
and sparse networks are primarily synergetic during early 
stages of development. However, in dense networks, they be- 
come primarily redundant by 30 DIV. Conversely, in sparse 
networks triplets are primarily synergetic by 30 DIV. Due to 
the fact that in vitro cultures have many similarities with in 
vivo tissue, we believe that the developmental features identi- 
fied here may also play a role in living organisms. These find- 
ings should be applicable to future research seeking to repli- 
cate or emulate brain functions and provide useful constraints 
for experiments of neural function. 
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